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We have searched for time-reversal symmetry breaking fields in the non-centrosymmetric super- 
conductor MgioIri9Bi6 via muon spin relaxation in zero applied field, and we measured the temper- 
ature dependence of the superfluid density by muon spin rotation in transverse field to investigate 
the superconducting pairing symmetry in two polycrystalline samples of signficantly different puri- 
ties. In the high purity sample, we detected no time-reversal symmetry breaking fields greater than 
0.05 G. The superfluid density was also found to be exponentially-flat as T— >0, and so can be fit to 
a single-gap BCS model. In contrast, the lower purity sample showed an increase in the zero-field 
/^SR relaxation rate below T c corresponding to a characteristic field strength of 0.6 G. While the 
temperature-dependence of the superfluid density was also found to be consistent with a single-gap 
BCS model, the magnitude as T— >-0 was found to be much lower for a given applied field than in 
the case of the high purity sample. These findings suggest that the dominant pairing symmetry 
in high quality MgioIrigBig samples corresponds to the spin-singlet channel, while sample quality 
drastically affects the superconducting properties of this system. 

PACS numbers: 76.75,+i, 74.70.Dd 



The spin-singlet Cooper pairs of conventional super- 
conductors require the corresponding superconducting 
order parameter to exhibit time-reversal symmetry [l(. 
For spin-triplet pairing, time-reversal symmetry is char- 
acteristic of only some of the possible superconducting 
states, but Anderson originally hypothesized that these 
should all have spatial inversion symmetry [2]. This then 
suggested that a material lacking an inversion center 
should be an unlikely candidate for spin-triplet pair- 
ing. However, a potential exception was found re- 
cently with the discovery of superconductivity in the non- 
centrosymmetric, heavy fermion compound CePt3Si[3] 
where the upper critical field was found to exceed the 
paramagnetic limiting field in this material, suggesting 
possible spin-triplet pairing. A mixed spin-singlet and 
triplet pairing state was proposed to resolve this contro- 
versy, and this exotic possibility has stimulated signifi- 
cant effort in both theoretical and experimental studies 
to understand superconductors lacking inversion centers. 

The proposed mixed pairing state was found to be 
a consequence of superconductivity in a material that 
exhibits spin-orbit coupling (SOC) but lacks inversion 
symmetry. In this case, the usually doubly-degenerate 
electronic bands become non-degenerate almost every- 
where in the Brillouin zone and the resulting super- 



conducting states can no longer be classified according 
to parity. While non-centrosymmetric, heavy fermion 
superconductors have been studied in more detail 
to-date, these materials often have magnetically-ordered 
states that compete or coexist with the superconduct- 
ing state, and such behaviour tends to significantly com- 
plicate penetration depth measurements. For example, 
the heavy fermion superconductor CePtaSi first becomes 
magnetically-ordered at 2.2 K before superconducting 
below 0.75 k(|, and the two states have been shown 
to coexist on a microscopic level Q- Transition-metal 
compounds such as Li 2 M 3 B (M = Pd, and Pt)[8|, |9|, 
MgioIrigBiedU, and M 2 Ga 9 (M = Rh, Ir)(llj are there- 
fore generally more straightforward systems for exploring 
the effect of SOC on superconductivity in the absence of 
inversion symmetry. 

The strength of the SOC in any material is determined 
by both the crystallographic structure and the elemental 
composition. The latter is well-illustrated through re- 
cent NMR[l3| and magnetic penetration depth[l3| mea- 
surements on Li2M3B, which found that increasing SOC 
by replacing Pd with Pt changes the superconducting or- 
der parameter from dominantly spin-singlet (I^PdaB) to 
nodal, spin-triplet (I^PtaB). Muon spin rotation mea- 
surements of the penetration depth in L^PdaB sup- 
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port the conclusion that L^PdaB is predominantly an 
s-wave BCS superconductor with only one isotropic en- 
ergy gap 

MgioIrigBi6 crystallizes in a non-centrosymmetric, 
body-centered cubic structure with space group I43m and 
lattice constant a=10.568 A[l5j]. A combination of re- 
sistivity, specific heat, and magnetization measurements 
indicate that this material is a superconductor with T c 
ranging between 4.5 and 5 K depending on the exact 
sample investigated ToL 16]. Since Ir is a heavy transi- 
tion metal SOC is expected to be strong in this material, 
and so questions about the superconducting pairing sym- 
metry are of significant interest. 

The superconducting pairing symmetry of 
MgioIrigBi6 has already been investigated through 
both heat capacity [la. Kl\ and photoemission 
spectroscopy [3, [li|, and the data is consistent with 
isotropic s-wave pairing in both cases. However, point 
contact spectroscopy[l6| and penetration depth mea- 
surements performed using a tunnel diode oscillator [20j 
suggest that the system is better described by two-gap 
superconductivity. The two-gap model is consistent with 
spin-split bands that are a direct consequence of the lack 
of inversion symmetry combined with significant SOC, 
and is expected in the case of a mixed pairing state with 
a significant spin-triplet component. To help resolve 
the apparent disagreement regarding the nature of the 
superconducting pairing symmetry in MgioIrigBi6, more 
experiments need to be performed. 

Muon spin rotation is a powerful technique for mea- 
suring the penetration depth and determining the na- 
ture of the pairing symmetry in superconductors. Due 
to maximal parity violation in the weak interaction de- 
cay of pions, the muons produced for such an experiment 
are initially ~100% polarized. Muon spin rotation (TF- 
^SR) employs a transverse-field geometry, as the muons 
are implanted into the sample one at a time with the 
initial polarization oriented perpendicular to the applied 
magnetic field. The muon spins exhibit Larmor preces- 
sion about the local magnetic field and then decay into a 
positron, plus two neutrinos which are not detected. The 
resulting positron travels preferentially in the direction of 
the muon spin at the time of decay, and so due to the 
asymmetric nature of this decay one can track the po- 
larization of the muon ensemble as a function of time by 
setting up positron counters around the sample. In type 
II superconductors, the mixed (or vortex) state gives rise 
to a spatial distribution of local magnetic fields; this field 
distribution manifests itself in the /iSR signal through a 
relaxation of the muon polarization. /iSR has been used 
successfully for some time to study the mixed state in 
superconductors; further details can be found in [ill, [jjj] . 

Zero field (ZF) /iSR is very complementary to TF-/iSR 
superconductor studies. The experimental set-up is es- 
sentially the same, but no field is applied. /xSR is an ex- 
tremely sensitive probe of magnetism (internal fields on 
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FIG. 1: Powder x-ray diffraction measurements of the two 
MgioIri9Bi6 samples used for the /J.SK measurements. The 
arrows indicate the dominant impurity peaks in the sample 1 
pattern. 



the order of 0.1 G can be readily detected), and so when 
investigating superconductors a ZF-/iSR experiment is 
generally performed in advance of TF-/1SR measurements 
to search for any anomalous magnetism that was not pre- 
viously detected by other techniques. This knowledge is 
quite important since the presence of magnetism in a 
superconductor can significantly complicate penetration 
depth measurements. 

ZF-/iSR is also useful to search for time-reversal 
symmetry-breaking fields; the presence of such fields lim- 
its the possible superconducting states and the associ- 
ated pairing symmetry. For example, time-reversal sym- 
metry is a prerequisite for any spin-singlet state, but 
this symmetry is broken for certain spin-triplet states, 
specifically for those which have a degenerate repre- 
sentation. The presence of two or more degenerate 
superconducting phases naturally leads to a spatially- 
inhomogeneous order parameter near the resulting do- 
main walls; this creates spontaneous supercurrents and 
hence magnetic fields near those regions. Another pos- 
sible origin of time-reversal symmetry-breaking fields is 
from intrinsic magnetic moments due to spin polariza- 
tion (for spin triplet pairing) and the relative angular 
momentum of the Cooper pairs . Broken time- reversal 
symmetry was previously detected in the unconventional 
superconductor Sr 2 Ru04 via /LtSR^24], as a local field of 
~ 0.5 G was found below T c at the muon site. This 
provided some of the strongest evidence that Sr 2 Ru04 is 
actually a p-wave superconductor. Similar time-reversal 
symmetry-breaking fields could be detected by fiSK if the 
spin-triplet contribution to the superconducting state in 
Mgi IrigBi 6 is significant. 

In this work, we have used both ZF and TF-/iSR to in- 
vestigate the nature of the superconductivity in two poly- 
crystalline Mgi IrigBi 6 samples of different quality. The 
former allowed us to search for time-reversal symmetry- 
breaking fields, and the latter allowed us to measure the 
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FIG. 2: ZFC and FC susceptibility measurements of the two 
MgioIrigBie samples used for the (iSR measurements. 



penetration depth and characterize the pairing symme- 
try. Our data from the high-quality sample is consis- 
tent with a single-gap s-wave model; this is in agree- 
ment with the earlier heat capacity [17] and photoemis- 
sion spectroscopy work (l8i fljj ]. 

MgioIrigBig samples were synthesized by a standard 
solid state reaction of pure Mg, Ir, and B elements 
as described elsewhere, and the chemical composition 
was accurately determined via electron diffraction (loj. 
Prior to the fiSK measurements, the samples were char- 
acterized by powder x-ray diffraction, performed on a 
Scintag XDS 2000 diffractometer with Cu-K Q radiation 
(A =0.15460 nm). Rietveld refinements of the structures 
were achieved using GSAS(25l. [2o|. and slightly different 
lattice constants of 10.5629(4) A and 10.5655(1) A were 
found for sample 1 and 2 respectively. This is likely due 
to more vacancies in the crystal structure of the former, 
and when one also considers that many more impurity 
peaks are observed in the corresponding x-ray pattern 
(shown in Fig. 1), this suggests that sample 2 is a much 
higher quality specimen with fewer impurities/defects. 
The magnetic susceptibility of each sample was also mea- 
sured with a SQUID magnetometer; the zero-field cooled 
(ZFC) and field-cooled (FC) measurements, shown in 
Fig. 2, confirmed nice superconducting properties. If we 
let T c correspond to the onset of the susceptibility drop, 
then T c ~ 4.8 K for sample 1 and 4.6 K for sample 2. 
The much sharper drop below T c , observed in both the 
FC and ZFC measurements of sample 2, provides further 
evidence for the better sample quality in this case. 

ZF-/iSR data was collected using the dilution refriger- 
ator in conjunction with the Ml 5 surface muon beam- 
line at TRIUMF. Fig. [3] depicts the ZF muon spin ex- 
ponential relaxation rate vs. temperature for our two 
MgioIrigBig samples. Note that the data was fit to the 



FIG. 3: ZF-fiSH relaxation rate vs. temperature. The relax- 
ation rate of sample 1 shows a small increase around T c , while 
the relaxation rate of sample 2 remains essentially constant 
within the error bars. Possible reasons for this discrepancy 
are outlined in the text. 

function 

A(t) = A e- xt G KT (t) (1) 

where Gkt{^) is the Gaussian Kubo-Toyabe relaxation 
function given by: 

G KT (t)= 1 - + 2 -(l-Ah*)e-^ 2t2 (2) 

The Gaussian Kubo-Toyabe relaxation is due to nuclear 
dipolar fields and is temperature-independent; we deter- 
mined this to be A = 0.207(4) /is -1 from a global fit 
incorporating all of our ZF data for the two samples. 
The exponential relaxation rate was set as a variable pa- 
rameter, and this is plotted as a function of temperature 
for the two samples in Fig. 3. Sample 1 shows a small 
increase in the exponential relaxation rate of ~0.05 (is 
below T c , corresponding to a characteristic field strength 
of B ~ 0.6 G. However, the exponential relaxation rate 
of sample 2 remains essentially constant below T c . Al- 
though we cannot completely rule out the presence of 
time-reversal symmetry breaking fields in sample 2, the 
current data does set an upper bound on the magnitude 
of these fields to be ~ 0.05 G. 

There are a few different scenarios that can explain 
the discrepancy between the ZF data of the two sam- 
ples. One possibility is that the relaxation increase 
in sample 1 is actually due to time-reversal symmetry- 
breaking fields. The presence of such fields suggests 
that the associated superconducting state is at least 
two-fold degenerate, and the spatially- varying order pa- 
rameter would then lead to the presence of supercur- 
rents in the bulk of the superconductor 23] . If this is 
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the case, one would generally expect these fields to ex- 
ist in sample 2 also, although we note that the do- 
main size may be substantially smaller in sample 1 due 
to the poorer quality. Many more muons would then 
stop near the domain walls in sample 1 and see a non- 
zero local field associated with the currents flowing at 
the domain boundaries; this could lead to the increased 
broadening below T c for this case only. Furthermore, 
NMR and electron diffraction measurements were per- 
formed very recently on two similar samples with com- 
positions Mg9.3iri 9 Bi 6 . 7 and Mg10.5Ir19B17.1p7]. The 
electron diffraction measurements clearly confirmed that 
the former was a better quality sample with fewer de- 
fects/impurities, and the NMR measurements indicated 
that only the other sample had a finite spin suscepti- 
bility in the zero temperature limit. These results then 
suggested that defects/impurities may enhance spin or- 
bit coupling in this system thus leading to an increased 
spin-triplet component of the superconductivity; this in- 
terpretation is also consistent with the observation of 
time-reversal symmetry-breaking fields in our lower qual- 
ity sample only. 

There are also other possible explanations for the dif- 
ference between the ZF data of the two samples not re- 
lated to time-reversal symmetry-breaking fields. For ex- 
ample, the increased ZF relaxation in sample 1 could be 
due to dilute magnetic impurities throughout the sam- 
ple. However, this seems rather unlikely since the onset 
of the increased relaxation rate occurs very close to T c 
and furthermore, these are not expected in this system. 
The other possibility is related to field-zeroing. If the 
applied field was not perfectly zeroed in our experiment, 
there could potentially be some trapped flux in the su- 
perconducting state of these samples. Sample 1 may be 
more efficient at trapping flux than sample 2 due to the 
lower quality of the former, and this could result in an 
increased relaxation below T c for sample 1 only. We note 
that the ZF measurements were conducted on both sam- 
ples back-to-back without ever applying a field during 
this time, and so any residual applied field should be the 
same for both samples. 

While it is not possible for us to distinguish between 
these different possibilities, we can conclude with cer- 
tainty that a high purity sample of MgioIrigBi6 showed 
no evidence for time reversal symmetry-breaking fields 
greater than 0.05 G. We also observed only very weak 
ZF relaxation in both samples at all temperatures, and so 
this allows for reliable penetration depth measurements 
to be performed. 

With this in mind we proceeded to collect TF-^iSR 
data, as measuring the penetration depth would allow us 
to learn more about the nature of the superconducting 
pairing symmetry. Note that all of the data was col- 
lected via field cooling. Fig. |4] shows the Fast-Fourier 
transforms (FFTs) obtained at both 6 and 0.1 K. Be- 
low T c , the characteristic lineshape of an Abrikosov vor- 
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FIG. 4: MgioIri9Bi6 sample 1 asymmetry vs. time and the 
corresponding fast-Fourier transforms at TF = 500 G for both 
6 K [shown in (a) and (b)] and 0.1 K [shown in (c) and (d)]. 
While the increased broadening is apparent for T = 0.1 K, the 
characteristic Abrikosov lineshape was not observed. This is 
likely due to the large background signal superimposed on the 
data. 

tex lattice was not resolved as a large background signal 
from the silver sample holder was superimposed on the 
data. However, we did find significant line-broadening 
below T c as expected for a superconductor (see Fig. 
2]). This broadening was analyzed as a Gaussian field- 
distribution; the width of this distribution a sc has been 
shown previously to be proportional to the superfluid 
density a sc oc n s /m* oc 1/A 2 [28]. a sc was first ob- 
tained by fitting the /iSR data in the time domain to the 
functional form: 

A(t) = A e _ff2 * 2 ' /2 cos(a;t + 0) + A bg cos(uj bg t + cj)) (3) 

The non-relaxing background term was included to ac- 
count for the muons stopping in the silver sample holder. 
The first term gives the total sample relaxation rate a; 
there are contributions from both the vortex lattice and 
nuclear dipole moments below T c . The contribution from 
the vortex lattice, a sc , was determined by quadratically- 
subtracting the background nuclear dipolar relaxation 
rate obtained from spectra measured above T c . 

Fig. [5] depicts the TF muon spin relaxation rate (pro- 
portional to superfluid density) vs. temperature for our 
two MgioIri9Bi6 samples. Sample 1 was investigated 
with both TF = 500 G and 3 kG, while sample 2 was 
investigated only in 500 G. The sample 1 data indi- 
cates that the superfluid density decreases with increas- 
ing applied field as expected. The superfluid density of 
the two samples at 500 G as T— > is noticeably dif- 
ferent; sample 2 has a higher superfluid density than 
sample 1 and this is consistent with the higher qual- 
ity of the former. However, the temperature depen- 
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FIG. 5: TF-piSR superconducting relaxation rate a sc vs. tem- 
perature. The data can be adequately fit to a single-gap BCS 
model, as indicated by the solid lines. 



dence of the superfiuid density is the same for both sam- 
ples: namely, the superfiuid density seems to saturate 
as T — > 0. Since conventional BCS theory predicts an 
exponentially-flat temperature-dependence of the super- 
fluid density at low temperatures, we proceeded to fit 
our data to this model. The agreement is quite good 
for both samples as seen in Fig. [5j and so this sug- 
gests that the pairing symmetry in this material can 
likely be described by a single-gap s-wave model. For 
TF = 500 G, fitting to this model gave us the follow- 
ing parameters: 2A/kT c = 4.74(9) (2A = 1.6 meV), 
T c = 3.8(1) K, and A(T— > 0) = 3570(20) A for sample 1 
and 2A/kT c = 4.93(5) (2A = 1.6 meV), T c = 3.7(1) K 
and A(T— > 0) = 3230(20) A for sample 2. Note that our 
2A/kT c values are only slightly larger than what one 
expects for weak-coupled BCS theory, and that the pen- 
etration depth values were obtained from a sc through a 
conversion procedure outlined in [29j . The significantly 
suppressed T c 's in only a modest applied field of 500 G 
are consistent with the relatively low H C 2 ~8 kG[l6[ of 
this system. 

Motivated by the earlier work proposing multigap su- 



perconductivity in MginIriqBifi|16L |20|. we also tried to 
fit our data to a phenomenological two-gap model as de- 
scribed elsewhere [13]. Since the fitting routine assigned a 
weighting of (within our error bars) to one of the gaps, 
we found that a two-gap model does not seem to describe 
our data unless the two gaps are very similar in magni- 
tude. However, we note that previous measurements of 
MgioIrigBi6 suggesting multigap superconductivity were 
essentially performed in ZF, while our /iSR measure- 
ments were performed in an applied field. While a re- 
cent /xSR penetration depth measurement performed by 



our group on the superconductor Ba(Feo.926Coo.o74)2As2 
found evidence for multigap superconductivity in that 
material, the same measurements also revealed that the 
relative weighting factor of the smaller gap quickly ap- 
proaches with increasing TF[3(|. This could also be 
applicable to the present case. 

Although the above results come with the warn- 
ing that previous fiSK measurements on polycrystalline 
cuprate samples have resulted in the correct temperature- 
dependence of the superfiuid density being masked, this 
may be less of a concern in the present case due to the cu- 
bic crystal symmetry and hence the isotropic penetration 
depth. In any event, future penetration depth measure- 
ments should be conducted on single crystals when they 
become available to confirm our results. 

In conclusion, we have used both ZF and TF-^SR to 
investigate the superconductivity in two samples of the 
non-centrosymmetric system MgioIrigBig. In the higher 
quality sample, ZF-^SR measurements find no evidence 
for time-reversal symmetry-breaking fields larger than 
0.05 G. TF-^iSR measurements of both samples yield a 
superfiuid density that is exponentially-flat at low tem- 
peratures, and so our data can be fit to a single-gap 
BCS model. These findings suggest that high quality 
MgioIi'igBi6 exhibits spin-singlet s-wave pairing, even 
though the lack of inversion symmetry and the large ex- 
pected SOC of MgioIrigBi 6 should yield a mixed pairing 
state, presumably with a large spin-triplet contribution. 
Our work also suggests that the superconducting prop- 
erties of this system are very sensitive to sample quality, 
and so further studies are required to investigate the roles 
of both SOC and defects/impurities on the superconduc- 
tivity of this material in more detail. 
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